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ABSTRACT 

We make use of deep HST, VLT, Spitzer and Chandra data on the Chandra Deep Field South to 
constrain the number of Compton thick AGN in this field. We show that sources with high 24/im to 
optical flux ratios and red colors form a distinct source population, and that their infrared luminosity 
is dominated by AGN emission. Analysis of the X-ray properties of these extreme sources shows that 
most of them (80±15%) are indeed likely to be highly obscured, Compton thick AGNs. The number 
of infrared selected, Compton thick AGNs with 5.8/im luminosity higher than lO'*^'^ erg s^^ turns out 
to be similar to that of X-ray selected, unobscured and moderately obscured AGNs with 2-10 keV 
luminosity higher than 10"*^ erg s~^ in the redshift bin 1.2-2.6. This "factor of 2" source population 
is exactly what it is needed to solve the discrepancies between model predictions and X-ray AGN 
selection. 

Subject headings: Active Galactic Nuclei 



^ . 1. INTRODUCTION 

Active Galactic Nuclei (AGN) are not only witnesses 
Q of the phases of galaxy formation and/or assembly, but 
are most likely among leading actors. Indeed, three semi- 
c/2 . nal discoveries indicate tight links and feedbacks between 
super-massive black holes (SMBH), nuclear activity and 
galaxy evolution. The first is the discovery of SMBH 
, in the center of most nearby bulge dominated galaxies, 
^ ■ and the tight correlation between their masses and galaxy 
' bulge properties (Gebhardt et al. 2000, Ferrarese & Mer- 
\^ , ritt 2000, Marconi & Hunt 2003 and references therein). 
OO ' The second is that the growth of SMBH is mostly due to 
accretion of matter during their active phases, and thcre- 
fore that most bulge galaxies passed a phase of strong 
Q nuclear activity (Soltan 1982, Marconi et al. 2004). The 
J — . third is that the evolution of AGN is luminosity depen- 
, dent, with lower luminosity AGN peaking at a redshift 
!• ' lower than luminous QSOs (Hasinger 2003, 2005, Fiore 
■ et al. 2003, Ueda et al. 2003, La Franca et al. 2005, 
K> , Brandt & Hasinger 2005, Bongiorno et al. 2007), a bi- 
5_j ■ modal behavior recalling the evolution of star-forming 
d ' galaxies and that of massive spheroids (Cowie et al. 1996, 
Franceschini et al. 1999, De Lucia et al. 2006). AU 
three discoveries imply that obtaining a complete cen- 
sus of accreting SMBH through the cosmic epochs and 
constraining accretion efficiency and feedbacks are cru- 
cial steps toward the understanding of Galaxy formation 
and evolution. 

First attempts to constrain models for the formation 
and evolution of structure in the Universe using the 
evolving optical and X-ray AGN luminosity functions 
have been presented by Granato et al. (2001, 2004), Di 
Matteo et al. (2005), and Menci et al. (2004, 2005). In 
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particular, the Menci et al. model links the evolution of 
the galaxies in the hierarchical clustering scenarios with 
the changing accretion rates of cold gas onto the cen- 
tral SMBH that powers the QSO (Cavaliere & Vittorini 
2000). The results of this model were encouraging, in the 
sense that it predicts a trend of lower luminosity AGN 
to peak at increasingly lower redshift, as observed. How- 
ever, from a quantitative point of view, the model over- 
predicts by a factor of about 2 the space density of low-to- 
intermediate luminosity (Seyfert like) AGNs at z=1.5-2.5 
with respect to present X-ray observations. Furthermore, 
Marconi et. al. (2004, 2007 in preparation) derived a 
SMBH mass function from the X-ray selected AGN lu- 
minosity functions (e.g. La Franca et al. 2005) that falls 
short by a factor of about 2 to the "relic" SMBH mass 
function, evaluated using the yisH — <^v I Mg^f— M^ 
relationships and the local bulge's luminosity function. 
The most likely explanation for both discrepancies is that 
present X-ray surveys, although very efficient to probe 
unobscured and moderately obscured AGN (with column 
densities up to a few lO'^"^ cm^^, the so-called Compton 
thin AGNs), miss most of the very highly obscured, but 
still strongly accreting objects, the so called Compton 
Thick AGNs, with a column density Njj > lO^'' cm~^ (see 
Comastri 2004). Indeed, only a handful of the faintest 
sources in the Chandra deep fields may be Compton thick 
(see La Franca et al. 2005 and Tozzi et al. 2006). So 
we still may be viewing just the tip of the iceberg of 
the highly obscured AGN population. Compton thick 
objects may well be more common at high redshift, as 
suggested on theoretical ground by Silk & Rees (1998) 
and Fabian (1999) and on observational ground by e.g. 
Gilh et al. (2001), Worsley et al. (2004, 2006), and La 
Franca et al. (2005). 

Compton thick AGN at z > 1 can be recovered thanks 
to the reprocessing of the AGN UV emission in the in- 
frared by selecting sources with AGN luminosities in the 
mid-infrared and faint near-infrared and optical emission 
(e.g. Martinez-Sansigre et al. 2005, 2006, Houck et al. 
2005, Weedman et al. 2006a, 2006b). We investigate fur- 
ther this issue making use of the multiwavelength data 
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Fig. 1. — \L\ spectral energy distributions of the four addi- 
tional templates of sources hosting highly obscured AGNs used in 
this work. From top to bottom: NGC6240; IRAS09104+41091; 
HELLAS2XMM A2690.75 and BPM16274_69 

obtained on the Chandra Deep Field South (Giacconi et 
al. 2002), one of the fields with the deepest coverage at 
optical, infrared and X-ray wavelengths. 

A i7o = 70 km s'^ Mpc-\ Om=0.3, Oa = 0.7 cosmol- 
ogy is adopted throughout. 

2. DATASETS AND SAMPLE SELECTION 

The selection and spectroscopic identification of com- 
plete AGN samples from mid-IR surveys is a rather dif- 
ficult task, because AGNs make only a small fraction of 
the full mid-IR source population. Using Spitzer IRAC 
colors results in samples significantly contaminated by 
star-forming galaxies (e.g. Lacy et al. 2004, Alonso- 
Herrero et al. 2006, Barmby et al. 2006, PoUetta et al. 
2006). Comparing the global observed Spectral Energy 
Distribution (SED) to AGN and galaxy templates proved 
to be more efficient in selecting AGN samples (Polletta 
et al. 2006, 2007). However, it is very difficult to assess 
the completeness of these samples given all complex se- 
lection effects. Furthermore, many obscured AGN may 
still be anyway missed by this technique. 

We adopt in this paper a somewhat different approach. 
We do not pretend to select all AGN through optical 
and infrared photometry. As explained before X-rays are 
much more efficient in selecting unobscured (i.e. broad 
line AGNs) and moderately obscured AGNs. We con- 
centrate our effort on highly obscured AGN only and 
we limit our analysis to the high infrared luminosity 
AGN population. The driving consideration is that dif- 
ferences between nuclear and star-formation emission are 
emphasized comparing the observed SEDs with galaxy 
templates over a range as broad as possible (also see 
Martinez-Sansigre et al. 2005,2006, Houck et al. 2005 
and Weedman et al. 2006, Magliocchctti et al. 2007). 
Our primary objective here is to validate our obscured 
AGN selection criteria and assess the magnitude of the 
corresponding selection effects using a careful analysis of 



the deep X-ray data available on the CDFS. 

2.1. The GOODS-MUSIC catalog 

We use in this paper the latest version of the GOODS- 
MUSIC catalog (Grazian et al. 2006). We hmit the anal- 
ysis to the region fully covered by deep VLT/ISAAC near 
infrared photometry (143.2 arcmin^) and to the sources 
with MIPS 24/^m fluxes F(24/xm)> 40/xJy (1729 sources). 
24/im fluxes for all objects in the catalog, were obtained 
following the procedures described in De Santis et al. 
(2007). 

The published version of the catalog contains z band 
and K band selected sources in the GOODS-south area. 
The revised version of the catalog used in this work in- 
cludes 46 objects that are detected only in the 4.5/im 
band, i.e. their z and K magnitudes are below the cho- 
sen detection threshold. Only 4 of these sources do not 
have counterparts in the optical and/or in the near in- 
frared images. At the flux limits adopted here, we do 
not detect any objects at 24/im that is not detected at 
shorter wavelengths. In the following we use for each 
entry of the catalog Vega magnitudes and cgs fluxes. 

The GOODS-MUSIC catalog includes both monochro- 
matic and total infrared luminosities, for the sources with 
a reliable spectroscopic or photometric redshifts. Total 
8-1000/im luminosities were computed by integrating the 
best fit galaxy and AGN templates. Monochromatic lu- 
minosities were computed by interpolating the observed 
SEDs at the rest frame wavelength of interest. A large 
(factor of 10-30) systematic uncertainty is associated to 
the total infrared luminosity, being this dominated by 
the contributions at wavelengths of ^ 100 — lOOO/im, 
well outside the infrared band used in this paper (1-24 
/im). Different models can produce similar fits below 
24/tm but give rise to large differences in the total in- 
frared luminosity. This systematic uncertainty is greatly 
reduced using the monochromatic luminosity at 5.8 /im 
(see e.g. Yan et al. 2007), since this wavelength is within 
the observed band up to z^3.2. For this reason in the fol- 
lowing we make use of the infrared luminosity at 5.8/im 
(ALA(5.8/im)) to characterize the infrared power of the 
CDFS sources. 

2.2. Galaxy and AGN templates 

A detailed fltting of the observed SEDs using AGN 
and galaxy templates was also performed. Templates in- 
clude passive galaxies, star-forming galaxies, unobscured 
AGNs and highly obscured AGNs. Spectral libraries in- 
clude both empirical templates (Coleman et al. 1980 
datasets, Polletta et al. 2007) and synthetic models 
(Bruzual & Chariot 2003, Fioc & Rocca-Volmerange 
1997). We put particular care on the description of 
highly obscured AGNs. In addition to the templates pre- 
sented by Polletta et al. (2007) we used four additional 
templates of sources hosting highly obscured AGNs (fig- 
urell]). The NGC6240 template includes U, B, V, J, H, K, 
IRAS and ISO photometry. The IRAS09104+41091 in- 
cludes SDSS u, g, r, i, z photometry, B, V, R, J, H, K pho- 
tometry, Spitzer IRAC and MIPS photometry, Spitzer 
IRS spectroscopy and IRAS photometry. For the two 
HELLAS2XMM sources we use the best fit SED model in 
Pozzi et al. (2007). These templates span a range of in- 
frared to optical flux ratio significantly broader than the 
obscured AGN templates used by Polletta el al. (2007). 
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Fig. 2. — Left panel: F(24/^m)/F(R) as a function of the 5.8/im luminosity for three X-ray source samples (GOODS-MUSIC, red symbols, 
ELAIS-Sl, green symbols and HELLAS2XMM, large blue symbols , Pozzi et al. 2007). Open circles = type 1 AGN; filled circles = non 
type 1 AGN; stars = photometric redshifts. Note that F(24/im)/F(R) of non broad line AGN is strongly correlated with the luminosity at 
5.8/im. Right panel: F(24;im)/F(R,) as a function of the X-ray to optical flux ratio X/O for the same source samples. All symbols as in 
left panel; skeleton triangles = sources without optical counterpart. 



2.3. Photometric redshifts and SED fittings 

Photometric redshifts were derived by Grazian et aL 
(2006) fitting only the part of the SED dominated by the 
integrated stellar population, i.e. A < 5.5/im (Grazian et 
al. 2006 and references therein). The synthetic models 
used are very accurate and complete in the treatment 
of the star-formation history and evolution of the stellar 
populations and of the dust content of the galaxy and its 
evolution. However, while dust extinction is easy to ac- 
count for because it is a line of sight effect, dust emission 
and reprocessing is much more complicated to model, 
since it depends largely on the assumed geometry and 
covering fraction. Dust emission was therefore not con- 
sidered in the models, and accordingly the bands above 
5.5/im rest frame are ignored in these fits (see Grazian 
et al. 2006 for further details). This approach produces 
robust photometric redshifts {Az/{l + z) < 0.05) for pas- 
sive galaxies, star-forming galaxies and obscured AGNs, 
where the nuclear optical and infrared emission is com- 
pletely blocked or strongly reduced by dust and gas along 
the line of sight. Unobscured AGNs with power law SEDs 
are excluded from this analysis (however, most of them 
have reliable spectroscopic rcdshift from Cimatti et al. 
2002, Szokoly et al. 2004, Le Fevre et al. 2004, Vanzella 
et al.2005,2006, Mignoh et al. 2005). 

To characterize each SED we repeated the fit to the 
observed SED with a library of empirical templates, fix- 
ing the redshift to the spectroscopic redshift or, if this 
is not present, to the photometric redshift obtained as 
described above. This library includes the 21 templates 
of passive galaxies, star-forming galaxies and AGNs by 
PoUetta et al. (2007) and the four templates of sys- 
tems hosting a highly obscured AGN in figure [T] The 
best fit template and normalization were used to mea- 
sure monochromatic and total infrared luminosities. 



2.4. Optical, near infrared and mid infrared color 
selection 

The longest and shorthest wavelengths at which deep 
photometry is available on the Chandra Deep Fields are 
the 24/im mid-infrared band covered by Spitzer MIPS 
and the optical bands covered by HST ACS. 24/im 
sources with faint optical counterparts must be either 
luminous AGN whose optical nuclear emission is blocked 
by dust and gas, or powerful dusty-starburst galaxies. 
The Mid-infrared to optical flux ratio (^) can therefore be 
considered a rough estimator of obscured activity (both 
nuclear and star-formation) in galaxies. 

Figure [2] shows F(24/im)/F(R) as a function of 
\Lx{5.8iJ,m), for three samples of X-ray sources. Un- 
obscured AGNs (open symbols in figure [2|) have 
F(24/L(m)/F(R) in the range 10-200, uncorrelated with 
XL\{5.8fiin), as expected because the nuclear emis- 
sion dominates both optical and mid infrared wave- 
lengths. Conversely, obscured AGNs (filled symbols) 
have F(24/im)/F(R) spanning a broader range, and fairly 
correlated with XL\{5.8fj,m). This behavior resembles 
that of moderately obscured AGNs detected in X-rays, 
for which the X-ray to optical flux ratio (X/O) is strongly 
correlated with the X-ray luminosity (Fiore et al. 2003, 
Eckart et al. 2006). The nuclear optical-UV light of 
these objects is completely blocked, or strongly reduced 
by dust extinction, and the optical fiux is dominated 
by the host galaxy. On the other hand, the 2-10 keV 
fiux is reduced by only small factors, even for obscur- 
ing gas column densities of the order of a few x 10^^ 
cm~^. X/O is therefore a good estimator of the ratio 
between the nuclear fiux and the host galaxy starlight 

^ F(24^im)/F(R); logF(R)=-0.4 X R- 22.5467. R magnitudes 
have been obtained by interpolating the V and I magnitudes in the 
GOODS-MUSIC catalog provided by HST/ACS 
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Fig. 3.— Left panel: F(24^jm)/F(R) as a function of the R-K color for two X-ray source samples (GOODS-MUSIC and HELLAS2XMM, 
large symbols, from Pozzi et al. 2007). Open circles = type 1 AGN; filled circles = non type 1 AGN; stars = photometric redshifts. 
Isodensity contours of all GOODS-MUSIC 24/xm sources are overlayed to the plot. The thick continuous lines mark the expectations 
of throe obscured AGN SEDs with redshift increasing from to 4 from left to right. The lower curve represents the colors of a typical 
low luminosity Seyfert 2 galaxy, the middle curve represents the colors of an obscured AGN from the Pozzi et al sample (A2690_75), 
the upper curve represents the colors of IRAS09104-I-41091. The black dashed lines are the expectations of the SEDs of two star-burst 
galaxies (M82, lower curve and Arp220, upper curve) for z=0-4. The dotted line is the expectation of a passive elliptical galaxy for z=0-4. 
Right panel: fraction of GOODS-MUSIC 24/xm sources as a function of the R-K color in three F(24^m)/F(R) bins: solid histogram = 
F(24^m)/F(R)< 200; dotted histogram = 200 <F(24^m)/F(R)< 1000; dashed histogram = F(24^m)/F(R)> 1000. 



flux. While the nuclear AGN X-ray luminosity can span 
several decades, the host galaxy R band luminosity has 
a moderate scatter, less than one decade, giving rise to 
the observed correlation between X/0 and the X-ray lu- 
minosity. Perola et al. (2004), Mignoh et al. (2004), 
Brusa et al. (2005) and Cocchia et al. (2007) found 
that the high X/0 sources tend also to be obscured in 
the X-rays, with column densities of the order of 10^^~^^ 
cm~^. It is therefore possible to conclude that a high 
X/0 ratio is a good indicator for both optical and mod- 
erate X-ray obscuration in high luminosity sources. In- 
terestingly, the F(24/im)/F(R) of X-ray selected sources 
is strongly correlated with X/O (figure [H right panel). 
This suggests that luminous Compton thick AGN, which 
are faint in X-rays because of the strong photoelectric 
absorption and Compton scattering, and cannot be se- 
lected using their X/0 flux ratio, can be recovered using 
the F(24^m)/F(R) ratio. 

Furthermore, since X-ray obscured AGNs tend to 
have red R-K colors (Brusa et al. 2005 and reference 
therein), one would also expect that Compton thick 
AGN have similarly red colors. Indeed figure [3] shows 
that the F(24^m)/F(R) of X-ray selected, obscured AGN 
is correlated with the R-K color, as expected. Fig- 
ure [3] also shows the iso-density contours of all the 
GOODS-MUSIC 24^m sources with F(24^m)> 40/xJy 
Intriguingly, the iso-density contours become narrow in 
F(24/im)/F(R) at high R-K values and extend toward 
the region occupied by obscured X-ray selected AGN at 
high F(24/um)/F(R) and high R-K values. The bimodal- 
ity at high values of F(24/im)/F(R) of the color distri- 
bution of the 24/Ltm selected sources is evident in the 
rightmost panel of figure [3l which shows the fraction of 



GOODS-MUSIC 24/im sources as a function of the R- 
K color in three F(24/im)/F(R) bins. While at low and 
intermediate F(24/im)/F(R) values the distributions are 
peaked at R-K~ 2.5 — 3.5 and decrease smoothly to- 
ward higher R-K values, the distribution of the sources 
with F(24/im)/F(R)> 1000 shows a strong excess at R- 
K> 4.5. 

It is interesting to note that most of the highly ob- 
scured AGN selected in the HELLAS 2XMM survey on 
the basis of their high X-ray to optical flux ratio (Pozzi 
et al. 2007) have F(24^m)/F(R) higher than a few hun- 
dred, and all have R-K> 4.5. Their SEDs are char- 
acterized by a passive early-type galaxy in the opti- 
cal and in the near infrared, and by an AGN compo- 
nent in the mid infrared. These SEDs redshifted up 
to z=4 are able to explain the extreme colors of the 
F(24^m)/F(R)> 1000 R-K> 4.5 sources, unlike the 
SED of even extreme star-forming galaxies like Arp220 
(see figure [3] left panel). This strongly suggests that 
most of the F(24/zm)/F(R)> 1000 R-K> 4.5 sources 
are powered by an active nucleus. Similar conclu- 
sions are found analyzing slightly different color dia- 
grams, like F(24^m)/F(R) vs. F(24/xm)/F(8/.tm) and 
F(24^m)/F(R) vs. F(3.6Aim)/F(z). 

3. X-RAY PROPERTIES OF EXTREME 24/iM SELECTED 
SOURCES 

Our candidate obscured AGNs are selected from the 
full 24^m GOODS-MUSIC sample with the request of 
having F(24/im)/F(R)> 1000 and R-K> 4.5. There are 
135 such sources. 

3.1. Sources with a direct X-ray detection 
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Fig. 4. — Left panel: the redshift distribution of the F(24/im)/F(R)> 1000 and R-K> 4.5 sources without direct X-ray detection (107 
sources) compared to those of the full GOODS-MUSIC 24/im sample (1649 sources) and of the X-ray selected AGNs in the GOODS-MUSIC 
area (150 sources). Right panel: the redshift-infrared luminosity {log{XL x{5.8fim))) plane for the full 24/xm GOODS-MUSIC source sample 
(isodensity contours); the CDFS X-ray AGN sample (open circles = type 1 AGN; filled circles = non type 1 AGN; stars = photometric 
redshifts); and the sample of F(24/im)/F(R)> 1000 and R-K> 4.5 sources without direct X-ray detection (filled triangles) 



Eighteen of the 135 sources have an X-ray detection 
in Alexander et al. (2003). Four other sources are not 
formally detected but have more than 4-5 counts (af- 
ter background subtraction) at the position of the 24/im 
source. In summary, 22 of the 135 sources (16%) have a 
significant X-ray emission directly visible in the Chandra 
images. 

Three of these 22 sources have a spectroscopic redshift, 
with narrow line optical spectra, the other 19 have a 
photometric redshift in the GOODS-MUSIC catalog (see 
Section 2); the median redshift and its interquartile range 
of these 22 sources are 2.1 (0.5; hereinafter interquartile 
ranges are indicated in brackets after the median values). 

The median monochromatic infrared luminosity at 
5.8/im is 44.42 (0.37). The X-ray luminosities are in all 
cases higher than 10^^ erg s^^, making them bona-fide 
AGNs. The median logarithmic ratio between the 5.8/L(m 
and the 2-10 keV luminosities is 1.07 (0.32). As a com- 
parison, the median of the same logarithmic ratio for the 
fuU GOODS-MUSIC X-ray sample (150 AGN with mea- 
sured redshift) is 0.69 (0.47). The probability that the 
two distributions are drawn from the same parent popu- 
lation is ^ 0.2%, using the Kolmogorov-Smirnov test. 

The hardness ratios indicate in most cases a hard, pos- 
sibly obscured. X-ray spectrum. Indeed, these sources 
are among the most obscured ones in the Tozzi et al 
(2006) analysis, all having column densities higher than 
a fewxlO^^ cm^^ and 2 having column densities as high 
as lO^'' cm~^. 

Table 1 gives the breakdown of the best fitting tem- 
plates (see Section 2) to the SEDs of these 22 sources. 
15 SEDs are best fitted by one of the templates in figure 

m 



3.2. Sources without a direct X-ray detection 



TABLE 1 

Table 1: template fits to the SEDs of the sources with 
F(24/im)/F(R)> 1000 and R-K> 4.5 



Template 


X-ray det. 


Not X-ray det. 


Ellipticals -|- SO 




2 


Spirals 




1 


M82-|-N6090-fArp220 




35 


119254 -1- Mark231 


4 


17 


Seyfert 1.8-2+red QSO 


1 


5 


A2690.75 + BPM16274_69 


5 


34 


IRAS09104+4109 


9 


11 


N6240 


1 


1 


Seyfert 1 +QSOs 


2 


1 


Total 


22 


107 



The total number of sources with F(24/im)/F(R)> 
1000 and R-K> 4.5 and no direct X-ray detection is 111 
(we excluded two sources which happen to lie within 5 
arcsec from an X-ray source). 

Four of these sources have a spectroscopic redshift, 99 
have photometric redshift in the GOODS-MUSIC cata- 
log. For four sources we could only compute a lower limit 
to the redshift. In conclusion, we have redshifts or limits 
for 107 sources. Both median redshift and infrared lumi- 
nosity are similar to those of the 22 sources with a direct 
X-ray detection. 

The redshift and infrared luminosities distributions of 
the F(24/.tm)/F(R)> 1000 and R-K> 4.5 sources are 
compared in figureHto those of the fuU GOODS-MUSIC 
24^m and X-ray selected samples. The F(24/im)/F(R)> 
1000 and R-K> 4.5 sources have a redshift distribution 
with median redshift 1.91 (0.30). Excluding the lower 
limits <z> = 1.9 (0.28). This distribution is shifted to- 
ward higher redshifts than both the full GOODS-MUSIC 
24/im source sample and the X-ray selected AGN sam- 
ple (figure d] left panel). The right panel of figure d] 
shows that moderately obscured. X-ray selected sources 
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Fig. 5. — Left figure: stacked X-ray images of 24/im GOODS-MUSIC sources. Images have sides of 18 arcsec, central circles have 2 arcsec 
radii. Upper panels refer to 111 sources with F(24/im)/F(R)> 1000 and R-K>4.5; lower panels refer to 73 sources with F(24/.4m)/F(R)< 200 
and R-K>4.5. Loft panels are stacks in the S 0.3-1.5 keV band, right panels are stacks in the H 1.5-4 keV band. Right figure: the hardness 
ratio (H-S) / (H+S) of the counts in circles of 2 arcsec radii as a function of the fraction of highly obscured AGN for the sample of 111 24;im 
GOODS-MUSIC sources with F(24/jm)/F{R)> 1000 and R-K>4.5. The solid curve is the result of Montecarlo simulations (see text for 
details); the two thick horizontal lines are the average hardness ratios measured for the F(24/xm)/F(R)> 1000 and R-K>4.5 (upper line) 
and F(24/im)/F(R)< 200 and R-K>4.5 sources (lower line). The colored bands mark the hardness ratio statistical uncertainties. 



are concentrated below z—1.5 and span a range of in- 
frared luminosities log(ALA(5.8/im)) ~ 42 — 45.3 Their 
median 2-10 keV and infrared luminosities are 43.16 
(0.56) and 43.84 (0.57) respectively. On the other hand, 
the F(24^m)/F(R)> 1000 and R-K> 4.5 sources are 
concentrated between z=1.2 and z=2.6 and have in- 
frared luminosities from log(ALA(5.8/im)) ^ 43.5 to 
\og{XLxi5.8^lm)) - 45.5, with a median 44.34 (0.30), 
luminosities similar to that of the X-ray selected AGNs 
at the same redshift. 

Table 1 gives the results of the template fits to the 
observed SEDs of the 107 sources with a redshift. In 
65% of the cases AGN templates provide the best fit. 
The majority of these SEDs are best fitted by one of 
the templates in figure [1] the others are best fitted by 
one of the AGN templates of PoUetta et al. (2007). 
Only one SED is best fitted by the template of a type 
1 AGN. 36 SEDs are best fitted by star-forming galaxy 
templates (in 30 cases the best fit is obtained using pow- 
erful star-forming galaxy templates, like those of Arp220 
and NGC6090). In conclusion, the results of the SED 
fitting analysis confirm that the majority of the sources 
with F(24^m)/F(R)> 1000 and R-K> 4.5 may be highly 
obscured AGNs. However, we remark that the results 
of the template fitting should not be always considered 
a quantitative determination of the nature of each sin- 
gle source, since in many cases different templates can 
produce fits with similar in the observed optical to 
24/im band, producing a degeneracy difficult to account 
for. Therefore, the results of the template fittings should 
only be taken as a qualitative indication of the popula- 
tion properties of the samples. 

3.3. X-ray stacking analysis 



To validate our highly obscured AGN selection sample 
and assess quantitatively the fraction of bona-fide AGNs 
in the F(24^m)/F(R)> 1000 and R-K> 4.5 sample, we 
performed a detailed "stacking" analysis of the X-ray 
data of the 24^m selected sources. Indeed, thanks to 
the low Chandra background, it is possible to increase 
the effective exposure time and derive average properties 
of undetected objects, using these stacking techniques: 
counts at the positions of known sources are co-added 
in order to probe X-ray emission substantially below the 
single source sensitivity limit. 

We performed a stacking analysis of the 111 sources 
not directly detected in the Chandra 1 Msec X-ray im- 
age. As control samples we used the 22 sources with 
F(24/Ltm)/F(R)> 1000 and R-K> 4.5 and an X-ray coun- 
terpart, a sample of sources with F(24^m)/F(R)> 1000 
and R-K< 3.5 (51 sources after the exclusion of the 
sources directly detected in the X-ray image) and a sam- 
ple of sources with F(24^m)/F(R)< 200 and R-K>4.5 
(73 sources after the exclusion of the sources directly de- 
tected in the X-ray image) . The total exposure times for 
the four source samples are 94 Msec, 19.5Msec, 43 Msec 
and 61 Msec respectively. 

The top panels of figure [5^) show the X-ray stack of 
the 111 24/xm sources with F(24/im)/F(R)> 1000 and 
R-K> 4.5 in two energy bands. This is compared to 
the stack of 73 sources with F(24^m)/F(R)< 200 and 
R-K> 4.5 and no direct X-ray detections (lower pan- 
els). We choose the bands 0.3-1.5 keV (soft band, S) 
and 1.5-4 keV (hard band, H) to keep the level of the 
internal background as low as possible and similar in the 
two bands. At a typical redshift of 2 these bands cor- 
respond to rest frame energies of 0.9-4.5 kcV and 4.5-12 
keV respectively. The stack of the high F(24^m)/F(R) 
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and high R-K sources produces a detection in both soft 
and hard bands, with count rates (0.98 ± 0.20) x 10^^ 
counts/s and (1.25 ± 0.26) x 10~^ counts/s in the two 
bands respectively (using a 2 arcsec radius extraction re- 
gion). The hardness ratio H-S/H+S measured for this 
sample is therefore 0.12±0.15. Conversely, no detection 
in either bands is obtained for the stack of 51 sources 
with high F(24^m)/F(R) and low R-K, with la upper 
limits of '--^ 2 X 10^^ counts/s. The stack of the low 
F(24/im)/F(R) and high R-K sources produces a signif- 
icant detection only in the soft band, with count rate 
(1.89 ±0.32) X 10^^ counts/s. The count rate measured 
in the hard band is significant at less than 3 a and corre- 
sponds to (4.9 ±1.9) X 10"^ counts/s. The correspond- 
ing hardness ratio is H-S/H+S=-0.58±0.18. The stack of 
the 22 sources with F(24/^m)/F(R)> 1000 and R-K> 4.5 
and an X-ray counterpart produces a hardness ratio H- 
S/H±S= -0.06±0.01. In conclusion, the stack of the 111 
sources with F(24/im)/F(R)> 1000 and R-K> 4.5 and 
without a direct X-ray detection produces a significant 
signal in both soft and hard X-ray bands. It is interest- 
ing to note that its hardness ratio suggests an average 
spectrum harder than even the average spectrum of the 
22 sources with similar IR and optical colors but with a 
direct X-ray detection. 

3.4. Simulations to assess the fraction of obscured AGN 
in the 24fJ-m source samples 

We used the observed flux in the stacked images, to- 
gether with the hardness ratio H-S/H+S, to constrain 
the fraction of highly obscured AGN in the samples. 
To this purpose, we generated simulated X-ray count 
rates and hardness ratios as a function of the fraction 
of the AGN, assuming that the F(24/im)/F(R)> 1000 
and R-K> 4.5 source sample is made by obscured AGNs 
and star-forming galaxies. We started from the observed 
redshift and infrared luminosities, and for the obscured 
AGN we assumed a log(ALA(5.8/xm)/L(2-10 keV)) lumi- 
nosity ratio chosen randomly in the range 0.4-1.2. The 
lower value is the typical ratio found by Silva et al. 
(2003) for low luminosity Seyfert 2 galaxies with col- 
umn density logNn S 24. It is also similar to the ratio 
found for the two HELLAS2XMM sources A2690_75 and 
BPM16274_69. The upper value is the ratio found for the 
powerful obscured QSO 1RAS09104+4109 (PiconceUi et 
al. 2007). 

For the star-forming galaxies we used a 
log(AiA(5.8^m)/L(2-10 keV)) luminosity ratio between 
2 and 2.8. These two values are obtained assuming a 
total infrared to 2-10 keV logarithmic luminosity ratio of 
3.6 (Ranalli et al. 2003) and correcting this for the ratio 
between the total infrared luminosity and the b.Sfim 
luminosity of powerful star-forming galaxies and spirals 
templates. 

We further assumed that the star-forming galaxies are 
not obscured in X-rays, while the AGNs are highly ob- 
scured. For the latter objects, we adopted a flat logN^f 
distribution from 23 to 26 cm~^. 

For a given fraction of AGN in the sample, we first 
decide, using a random generator, whether each source 
is an AGN or a star-forming galaxy. Then, if the source 
turned out to be an AGN, we choose an absorbing column 
density using again a random generator and the assumed 
logNjj distribution. We computed unobscured 2-10 keV 



luminosities from infrared luminosities, and then X-ray 
fluxes by folding a power law spectrum with spectral en- 
ergy index of 0.8 reduced at low energy by photoelectric 
absorption with a Chandra response matrix. For column 
densities ^ 10^^ cm^^ we assumed that the direct emis- 
sion is completely blocked by photoelectric absorption 
and Compton scattering. For these sources we assumed 
a reflection component with normalization 1/100 of the 
direct component and same spectral index. 

A small fraction of the simulated count rates (20-30%) 
is higher than the Chandra detection limit. This was ex- 
pected, given the fact that 22 out of 135 sources with 
F(24^m)/F(R)> 1000 and R-K> 4.5 do have X-ray 
counterparts. Since we are interested in the fraction of 
AGN in the sources not directly detected in the Chandra 
images, these count rates have been excluded from the 
following analysis. 

The result of the simulations is shown in figure [5l right 
panel. According to our simulations, the average hard- 
ness ratio of the F(24/im)/F(R)> 1000 and R-K> 4.5 
source sample without X-ray detection is reproduced if 
80±15% of the sources in the sample are highly ob- 
scured AGNs. The hardness ratio of the sources with 
F(24^m)/F(R)< 200 and R-K> 4.5 is reproduced if the 
fraction of the obscured AGN is in the range 0-20%. 
Changing the assumptions made to produce the simula- 
tions within reasonable ranges changes only slightly this 
result. 

The result of the simulations concerning the 
F(24^m)/F(R)> 1000 and R-K> 4.5 source sample is 
well consistent with the indication coming from the 1- 
24/im SED fitting with galaxy and AGN templates pre- 
sented in the previous section. 

4. DISCUSSION AND CONCLUSIONS 

We selected a source sample with extreme mid-infrared 
to optical flux ratio (F(24/im)/F(R)>1000) and red opti- 
cal colors (R-K>4.5) in the GOODS CDFS area. These 
sources are among the most luminous sources in the 
GOODS 24/im sample, because of the correlation of 
F(24/im)/F(R) with the infrared luminosity (see figures 
[2] and 31 right panel) . The fraction of these sources not 
directly detected in the Chandra images produces a sig- 
nificant stacked signal in both Chandra soft and hard 
X-rays bands. A detailed analysis based on Montecarlo 
simulation shows that the stacked count rates and hard- 
ness ratios can be reproduced if this source population is 
dominated by highly obscured (Nh > a few x 10^^ cm~^) 
AGNs. 

This conclusion is further confirmed by the following 
consideration. If we assume that this source popula- 
tion is dominated by star-forming galaxies, we can use 
the observed infrared luminosity and the dust-corrected 
UV luminosities to derive two different estimates of star- 
formation rates. On one hand total infrared luminosi- 
ties have been converted to a star-formation rate us- 
ing the following formula: SFR = Lmergs^^ /2.2A x 
10"^ Mo/yr, (Kennicutt 1998). On the other hand the 
UV star-formation rate is an output of the fits of the 
observed SEDs with synthetic models. These models 
have been obtained using the Bruzual & Chariot (2003) 
models, parametrizing the star-formation histories with 
a variety of exponentially declining laws (with timescales 
ranging from 0.1 to 15 Gyrs) and metallicities (from Z = 
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O.O2Z0 to Z = 2.5Zq, see Grazian et al. 2006 and refer- 
ences therein for further detail). Dust extinction is added 
to all star-forming models, with < E{B — V) < 1.1 and 
a Calzetti attenuation curve. Given the very red colors 
of our objects in the optical bands, they are typically fit 
with star-forming templates with a large dust extinction 
{E{B — V) > 0.5 ), resulting in a star-formation rate 
much larger than what derived by a simple, dust free 
UV-to-SFR conversion. 

The median logarithmic star-formation rate from in- 
frared luminosity is 2.31 (0.59) Mo/yr. The median log- 
arithmic UV star-formation rate is 0.83 (O.8O)M0/yr, a 
factor 30 lower. The large mismatch between these two 
estimates strongly suggests that the infrared luminosity 
of this source sample is not dominated by star-formation 
but rather by accretion. As a comparison, the median 
infrared and UV star-formation rates of the sources with 
F(2Vm)/F(R)< 200 and R-K>4.5 in our control sam- 
ple are 1.26 (0.38) and 1.10 (0.30), consistent with each 
other, suggesting that the infrared luminosity of these 
sources is not dominated by nuclear accretion (in agree- 
ment with their X-ray hardness ratio, see previous sec- 
tion). 

We can now compare the number of highly obscured 
AGNs selected at 24^m and not detected in X-rays, 
with the number of unobscured and moderately obscured 
AGNs directly seen in the X-ray images. To reduce the 
importance of complex selection effects we limit this anal- 
ysis to the sources in the redshift bin 1.2-2.6. At z=2.6 
the 5.8/im luminosity corresponding to the limit of 40/iJy 
at 24^m is log(AiA(5.8^m))=44.3-44.2 for Seyfert 2,1.8 
galaxy SEDs, and 44.4-44.5 for the four SEDs in figure 
[TJ Assuming a log(ALA(5.8/im)/L(2-10 keV)) luminosity 
ratio in the range 0.4-1.2 implies unobscured 2-10 keV lu- 
minosities in the range 43.3-43.9. This suggests to limit 
the comparison to the sources with logL(2-10 keV)> 43 
and log(ALA(5.8Atm)) > 44.2. 

In the 1.2-2.6 redshift bin there are at least 
46±10 sources with F(24/.tm)/F(R)>1000, R-K> 4.5, 
log(ALA (5.8/im)) > 44.2 and no direct X-ray detection 
(80±15% of 57 sources), which are probably highly ob- 
scured AGNs, and 44 X-ray selected AGNs, 7 of which 
with broad lines in their optical spectra. The number 
of 24/im selected, presumably highly obscured AGNs, 
missed by the CDFS X-ray survey is therefore similar 
to the number of X-ray selected AGNs. 

The median 5.8/im luminosity of the 57 
F(24/im)/F(R)>1000 and R-K> 4.5 sources with- 
out direct X-ray detection in the redshift and luminosity 
bins defined above is log(ALA(5.8/im))=44.48 (0.30), 
implying a median unobscured 2-10 keV luminosity in 
the range 43.3-44.1, adopting our assumption on the 
infrared to X-ray flux ratio. At z=2.6 these luminosities 
translate to X-ray fluxes between a few times 10~^^ 
and 10~^^ erg cm~^ , meaning that these sources 
would have been easily detected if they were not highly 
obscured. 

We have compared our findings with the predictions of 
the La Franca et al. (2005) and GiUi et al. (2007) models. 
In the redshift bin 1.2-2.6 and in the area covered by 
our 24/im sample the La Franca et al. (2005) luminosity 
function predicts about 75 AGNs with logL(2-10 keV)> 
43, 20 of which with a column density ^ lO^'' cm^^. The 



absorption distribution adopted by the Gilli et al. (2007) 
AGN synthesis model for the Cosmic X-ray Background 
predicts about 40 Compton Thick AGNs with the same 
limits for redshift and luminosity. We find 46±10 24/im 
selected, highly obscured AGNs, a number about twice 
the La Franca et al. (2005) prediction but well consistent 
with the GilH et al. (2007) prediction. 

Marconi et. al. (2004, 2007 in preparation) derived 
a SMBH mass function from the X-ray selected AGN 
luminosity functions that falls short by a factor of about 
2 to the "relic" SMBH mass function. This difference can 
be greatly alleviated adding the population of infrared 
selected, highly obscured AGNs to the AGN selected in 
the X-ray band below 10 keV. 

Finally, we note that the population of infrared se- 
lected, highly obscured AGNs, can help in reconciling 
the predictions of models of galaxy evolution with the 
observed AGN number densities at z^l — 1.5. Indeed, 
the Menci et al (2004) model predicts a number of low-to- 
intermediate luminosity AGNs at z=1.5-2.5 about twice 
with respect to that measured through X-ray 2-10 keV 
surveys (La Franca et al. 2005, Fiore 2006). We found 
above that in the redshift bin 1.2-2.6 the number of X- 
ray selected AGN with logL(2-10keV)> 43 is comparable 
to that of the infrared selected, highly obscured AGNs. 
Therefore, adding together the two populations would 
give a total number density of AGNs similar to that pre- 
dicted by the Menci et al. (2004) model. More detailed, 
quantitative, AGN number density calculations, span- 
ning wider redshift and luminosity ranges, are beyond 
the purpose of this paper and will be presented elsewhere 
(Fiore et al. in preparation). 

Most of the GOODS-MUSIC 24/im selected sources 
with high F(24/im)/F(R), high R-K and no X-ray detec- 
tion have a very faint optical and near infrared counter- 
part, only 14% of the sample having R< 26. Furthermore, 
only a handful of objects have a mid infrared flux high 
enough (F(24/imj; 0.5mJy) to allow Spitzer IRS spec- 
troscopy. The spectroscopic identiflcation of ^^90% of 
the sample must therefore await the advent of ELTs or 
the launch of JWST. 

Our analysis is limited to AGNs of intermediate lu- 
minosity at z~ 2. To extend the coverage of the 
luminosity-redshift plane requires a complementary ob- 
servation strategy. Particularly useful to this purpose is 
the SWIRE survey, which has both Spitzer and optical 
medium-deep coverage on 50 deg^ of sky. As an ex- 
ample, the SWIRE survey contains hundreds of sources 
with extreme F(24/tm)/F(R) flux ratios and red optical- 
near infrared colors. Interestingly, several dozens of these 
extreme sources have 24/im flux higher than ~ 1 mJy, al- 
lowing Spitzer IRS spectroscopy, and/or optical magni- 
tude brighter than R'^ 25, allowing spectroscopy with 8m 
class telescopes. Furthermore, all these sources are well 
within the reach of Herschel instruments between 75 and 
500 /im. Such long wavelengths observations can greatly 
help in separating nuclear activity from star formation 
when the two components have comparable integrated 
luminosities, e.g. in low luminosity, Seyfert like, AGNs. 
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